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Abstract 

We consider an energy harvesting transmitter broadcasting data to two receivers. Energy and data arrivals are 
assumed to occur at arbitrary but known instants. The goal is to minimize the total transmission time of the packets 
arriving within a certain time window, using the energy that becomes available during this time. An achievable 
rate region with structural properties satisfied by the two-user AWGN BC capacity region is assumed. Structural 
properties of power and rate allocation in an optimal pohcy are estabUshed, as well as the uniqueness of the 
optimal policy under the condition that all the data of the "weaker" user are available at the beginning. An iterative 
algorithm, DuOpt, based on block coordinate descent that achieves the same structural properties as the optimal 
is described. Investigating the ways to have the optimal schedule of two consecutive epochs in terms of energy 
efficiency and minimum transmission duration, it has been shown that DuOpt achieves best performance under the 
same special condition of uniqueness. 

Index Terms 

Packet scheduling, energy harvesting, AWGN broadcast channel, energy-efficient scheduling. 
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I. Introduction 

The basic offline problem of energy-efficient packet transmission scheduling [2], [3], [4], [5] is to 
assign code rates (consequently transmission durations) to a set of packets whose arrival times are known 
beforehand, so that they are all transmitted within a given time window with minimum total energy. The 
solution needs to strike a tradeoff between energy and delay based on the observation that energy per bit 
with many ideal and suboptimal coding schemes is convex and monotone increasing with rate. Recently, 
the problem has been reformulated with a model where energy gets "harvested" or replenished at certain 
known instants [6]. 

While in the former formulations transmission rate needs to be adapted to the arrival rate of information 
here it is adapted to the generation rate of energy. Considering both of theses effects introduces a richness 
to the problem on top of the initial model. The point-to-point problem in [6] was recast for finite energy 
storage [7] and for a wireless fading channel [8]. The formulation has been extended to an AWGN BC 
in [9], [10], considering a static pool of data to be sent at the beginning of the schedule. The same BC 
problem was also studied under a limited battery constraint [11]. 

The problem in [9] and [10] is reformulated in [14] relaxing the assumption that data is ready at the 
beginning of the schedule. This paper, extending the work in [14], considers the broadcast problem where, 
given an average transmit power constraint, rates are picked from an achievable rate region which obeys 
certain structural properties satisfied by the AWGN BC. The sender (transmitter) gets replenished with 
arbitrary amounts of energy as well as data packets of arbitrary length destined to each user at arbitrary 
points in time. 

The choices of power level and the rates to individual receivers across time is called a schedule. An 
optimal scheduling policy is defined to be one that transmits all the bits that have arrived within a certain 
time window, in the minimum possible amount of time 7"°?*. The policy is allowed to use as many energy 
harvests as it needs, provided it respects causality (no energy is used before it is harvested.) The problem 
considered in this paper is an offline problem, where data arrival and energy harvest instants and amounts 
are assumed to be known in advance. Although this kind of prior information of data and energy arrivals 
may not be a widely applicable assumption to real-world problems, results obtained from this work help us 
understand the nature of an optimal solution and boundaries on the best performance. Online formulations 
have also appeared in the literature. Notably, [12], [13] develop online scheduling policies for multihop 
networks on finite-horizon and infinite horizon problem formulations, respectively. 
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To minimize the overall transmission duration, rates to the individual users should be chosen in a way 
such that transmission is fast and energy efficiency is satisfied. However, as we choose higher rates, we 
loose from energy efficiency. Balancing between fast and energy efficient transmission, the decision of 
rates needs to be based on the sequence of energy harvests and data arrivals. 

This paper essentially bridges the work that considered scheduling on a Broadcast Channel data that 
becomes available at arbitrary points in time [2] and work that considered chunks of energy becoming 
available at arbitrary points in time [9]. It can also be viewed as the multiuser correspondent of the second 
problem considered in [6]. The main contributions are: (1) to establish structural properties of the optimal 
schedule, (2) the uniqueness of the optimal policy when all weaker user data is ready at the beginning, 
(3) an iterative algorithm which returns the optimal schedule under the same condition. It is shown in 
Section n that in an optimal policy, transmit power used is constant within each epoch, and may only rise 
from one epoch to the next, so that once it starts, the transmitter never lowers its power until it finally 
goes silent. On the other hand, the transmitter should increase its power only under certain conditions. 
These conditions, along with other structural properties of power and rate are established in Section III. 
Next, the uniqueness of the optimal policy is established under the condition that all of the weaker user's 
bits are available at the beginning. Finally, an iterative algorithm (that we refer to as DuOpt) based on the 
nonlinear block descent method which returns a feasible schedule carrying the same structural properties 
that the optimal is shown to have, is described. It has been shown that DuOpt returns the optimal schedule 
in case of static pool of weaker user data at the beginning of the schedule. We start by giving the problem 
statement in the next section. 

II. System Model 

Consider a broadcast channel with one transmitter and two receivers. Arbitrary amounts of energy, 
{Ei < oo, i = 1, 2, . . .}, as well as data for each user {Bf\Bf''' < oo, i = 1, 2, . . .} become available to 
the sender at arbitrary times U. A possible sequence of data and energy arrivals is illustrated in Fig. 1. E{t) 
denotes the total energy that has been harvested in [0,t) (regardless of how much of it has been used.) 
Similarly, Bi{t) and B2{t) denote the total number of bits destined to the first and second user, respectively, 
that the sender has obtained in [0, i). The interval between any two sequential arrival events (regardless 
of energy or data) will be called an inter-arrival epoch. The length of the i*^ epoch is = tj — tj_i. 

In this offline problem, all the future arrival times and amounts of energy and bits are known by the 
sender at i = 0. It is also assumed that harvested energy and data are available for use instantaneously 



as they arrive, and code rate and transmission power decisions can be changed instantaneously. However, 
codeword block lengths will be chosen such that each codeword is sent completely within a single epoch 
(note that starting and ending times of epochs are known ahead of time), so that no arrival event occurs 
during a codeword. Consequently, the power and rate pair decision will be fixed throughout each codeword. 
We are interested in minimizing the total transmission time for packets arriving by a certain time 
< oo, so W.L.O.G., set Bi{t) = Bi(W) for t > W, i = 1, 2. A schedule, which is a sequence of power 
and rate allocations, is feasible if it sends Bi{W) < oo bits to the 1** user and B2{W) < oo to the 2"*^ 
user (with a certain level of reliability^), without violating causality (at any time, using available energy 
and data by that time). We are interested in finding among all feasible schedules one with the smallest 
completion time, T°p*. 

The structure of the achievable rate region will be based on the two-user AWGN BC. The capacity 
region of a two-user discrete time AWGN BC with average power constraint P, noise variance u^, where 
the user's channel gain (si > 0) is larger than the 2"'^ user's (s2 > 0), consists of rate pairs (ri,r2) 
satisfying: 

where a, (0 < a < 1), denotes the fraction of P used for the 1** user. Since si > S2, the I''* user will 
be referred as the "stronger user", and the 2"*^ as the "weaker user". From (1) each user's rate can be 

expressed as a function of the other's and P, as ri = hi(P,r2), r2 = h2{P,ri). The rate functions hi and 
/i2 defined on 3^+ x 5^+ will be assumed to satisfy the following properties: 

1) Nonnegativity: hi{P,r) > 0,h2{P,r) > 0. 

2) Monotonicity: hi{P,r), h2{P,r) are both monotone decreasing in r, and monotone increasing in P. 

3) Concavity: hi{P,r), h2{P,r) are concave in P and r: ^'gff''"^ < 0, ^'^^'"^ < 0, for i e 1,2. 

d^hi{P,r) > p) d^hi(P,r) ^ 
drdP — dPdr — 

d^h2{P,r) _ r. dVi2{P.r) _ A 
•^^ drdP dPdr ~~ ^■ 

The results in the rest of the paper will be valid for any rate function satisfying (l)-(5), which are also 
satisfied by the AWGN BC [9]. 

It is straightforward to show that one can restrict attention to feasible schedules that do not change 
their power and rate allocations within epochs. 

'The achievable rate regions will be impUcitly assumed to correspond to a certain constant tolerable error probability respecting which it 
is possible to transmit a finite number of bits with a finite amount of energy per bit. 



Lemma i; In an optimal schedule, the power and rate pair remain constant within all epochs, except 
for the epoch during which the schedule ends. 

Proof. During an epoch, there are no energy or data arrivals and the claim is identical with the one stated 
and proved in Lemma 2 of [9]. The power will drop to zero when the schedule ends, which is in general 
within (and not necessarily at the end of) the last epoch used by the schedule. ■ 
With this, we will take rate and power assignments constant during an epoch. Let Pj be the total 
transmit power and rji be the rate assigned to the j*'* user during epoch i. Similarly, Pji represent the 
power assigned to user during epoch i. We are now ready to state the problem in terms of power 
and rate allocations to epochs, more precisely, an assignment of power and the stronger user's rate to 
each epoch (the weaker user's rate is thus determined). A final technical assumption will be useful 
in stating the problem: we shall assume that there is some k^"^ < oo such that there is at least one 
feasible schedule that ends within the first epochs. In other words, k^^ is an upper bound for epochs 
to be considered. In problem statement, k* denotes the last epoch of an optimal schedule, where k* < /c'^p. 

Problem 1: Transmission Time Minimization of Data Arriving at Arbitrary Points on an Energy 
Harvesting BC: 

Minimize: T = T{{Pi,ru}i<i<k^v) 

subject to: Pi > , < m < /ii(Pi,0) , r2i = hiiPuVu) 

k 

Y,Piii<E{tk) 

k* k' 

Pi^i + Pik'+i) {T-Y.^i)< E{T) (2) 

i=l i=l 

k k 

^rnC, < Si(ife) , Y.r2^i^ < B2{tk) (3) 

i=l i=l 

i 

for k = 1,2, k* = max{z : < T} 

fc* fc* 

i=l i=l 

k* k" 

^r2.e^+r2(fe.+l)(T-^Ci) = B2(T) (4) 

i=l i=l 

We will refer to (2) and (3) as energy and data causality constraints, respectively, as these ensure 
no energy is consumed and no bit is transmitted before becoming available. In addition, when the A;*'* 
inequality in (2) holds with equality, we shall say that A;*'* energy constraint is active. Similarly, equality 
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case in (3) will be referred as a data constraint being active. Finally, the feasibility constraint (4), ensures 
all the data bits destined to each user are transmitted. 

In the next section, we investigate structural properties that any optimal schedule has to satisfy. 

III. Structure of an Optimal Policy 

Lemma 1 recorded that in an optimal schedule power can only change upon a data arrival or energy 
harvest. The next result states that when power changes, it can only increase. The key to the proof is that 
more "bits per joule" can be sent by evenly distributing energy across a time interval (i.e., maintaining a 
constant power level, which is a consequence of the convexity properties of our rate functions.) If an even 
distribution of power requires transferring energy or bits to the latter epoch, it can always be done; hence, 
total transmit power never decreases in time. But, power may increase in time, because even distribution 
of power may result in unmet causality constraints. We state these results in Lemma 2. 

Due to space constraints, the proofs of the following results (Lemma 2 through Lemma 4) are omitted, 
and given in [18]. 

Lemma 2: (For proof see [18]) Consider an optimal schedule that ends during epoch k*. Power is 
non-decreasing with epoch index, i.e. Pi < Pj+i for i — 1,2, ... k* — 1. 

As stated in Lemma 2 power cannot decrease, yet may rise in time. In the next Lemma we note what 
is necessary condition for such a rise to occur in an optimal policy. 

Lemma 3: (For proof see [18]) In an optimal policy, power can only rise at ti (end of epoch i) if at 
least one of the conditions below holds: 

a) Energy constraint is active at point ti. (i.e.the ith energy constraint is active) 

b) The data constraints for both users are active at point ti. (i.e., the set of constraints in (3)) 

c) The weaker user's data constraint is active and data arrival to the weaker user occurs at time ti. 
The next set of results illustrate the structure of rate allocation in conjunction with the power allocation 

in an optimal policy. 

Corollary 1: (For proof see [18]) In an optimal policy, 

1) If power increases upon a data arrival for the second user, data to be sent to the weaker user have 
been finished by this event. 

2) If power rises upon a data arrival for the stronger user, all available bits have been sent by this event. 

3) If power increases upon an energy harvest, all energy available at the beginning of the former constant 
power band has been consumed by this energy harvest. 
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In the rest, some properties will be proved under the condition that all weaker user data is available at 
the beginning. We shall abbreviate this condition as follows: 

Definition 1: Weaker User Full Buffer Condition (WUFBC) is said to be satisfied whenever all of the 
data of the weaker user is available at the beginning of transmission. That is, B2{W) = i?2(0). 

The following lemma states an important feature of the stronger user rate distribution under WUFBC. 

Lemma 4: (For proof see [18]) Consider two consecutive epochs i and i + 1 of a given schedule, 
ending at U and U+i by definition, and suppose WUFBC holds for the problem instance. The following 
is necessary for the rate and power allocation to these two epochs of the given schedule to be locally 
optimal: The stronger user's rate is constant throughout [ti-i,ti), and Furthermore, the rate may 

jump up at i = U (staying constant otherwise) if at least one of the below is true: 

1) There is data arrival to the stronger user at t = and all the data that arrived before t = ti has been 
transmitted by U. 

2) An energy harvest occurs at t — ti and all of the power has been used for the stronger user during 
epoch i. 

We investigate the unique solution of Problem 1 in the next section. 

IV. Uniqueness of the Optimum Schedule Under WUFBC 

In the following lemma we note that an optimal schedule uses all energy harvested by the time the 
schedule ends completely. 

Lemma 5: The energy consumed by an optimal schedule that ends at T°p* is equal to £'(T°p*). 
Proof. To reach contradiction, consider an optimal schedule that consumes less energy than it harvested 
and has leftover energy in its energy buffer at 7"°^*. The remaining energy in the buffer could have been 
used in the last epoch to decrease the transmission completion time, which contradicts the minimality of 
T°P*. Hence, this schedule cannot be optimal. ■ 

Next, we show the uniqueness of the optimal schedule under WUFBC. 

Theorem 1: There is a unique optimum schedule under WUFBC, i.e., a unique power-rate allocation 
achieving T°p*. 

Proof. Suppose that there are two distinct optimal schedules, and S^, which have equal power and 
rate assignments until t^-i and differ for the first time at epoch s. Consider that the corresponding 
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power allocation vectors, P"^ and P^, also differ at epoch s such that Pf^ — Pf^i e {l,2,..,s — 1} 
and Pf^ < P^. First, assume that P"^ remains constant after epoch s, i.e., Pf = P^^i > s. By 
definition, both schedules end at T°p*. The total energy consumption of would be less than that of 

by r°P\ i.e., (E?:i^^^ei + (7^°P*-ife*+i)n^H-i) < {EtlP^^C^ + (T^"' - tk*+i)P,^*+,) , which 
contradicts Lemma 5. Hence, total transmit power of 5"^ cannot remain constant after tg. Since total 
transmit power is nondecreasing (See Lemma 2), it should increase after epoch s and before the end of 
transmission, i.e., P^ < , 3m e {s, s + 1, Since there are no data arrivals for the weaker 

user, the increase in total transmit power is either due to energy constraint being met or due to all 
the packets arrived by the time tu having been transmitted (cf. conditions (a) or (c) in Lemma 3). As 
Yl^=s ^t^i < Sr=s ^Fii'' "^^^ consumed all the available energy at the end of epoch u. Hence, 

must have transmitted all the bits arrived until which means that has transmitted at least the same 
number of bits to both users while consuming less energy than between and tu, which contradicts 
the optimality of . Therefore, if there are two distinct optimal schedules, S"^ and , their power 
allocation vectors cannot be different, i.e., P^ = P^. 

Now, consider two rate pair vectors, R"^ and R^, where {r^^, r^) = (r^, r^),Vie{l,2,..,s — 1} and 
r^^ < fis- Let the rate of the stronger user in 5*^, {rj^} stay constant after ts-i- By Lemma 4 rate of 
the stronger user cannot decrease, hence the rate of the stronger user in would be larger than that 
of S"^ after epoch s, i.e., = < rf^ < rfj , Vj e {s, s + 1, ...,k — 1}. Since both schedules 

end transmission at the same time, transmits fewer bits to the stronger user than does, which 
contradicts the fact that optimal schedule transmits all the packet arrivals by the end of transmission. 
Therefore, the rate of the stronger user in S"^ cannot stay constant after epoch s. Now suppose that rate 
of the stronger user in S"^ increases at the end of epoch u, i.e., r^^ < r^^+i) , 3it e {s, s + 1, ...,k*}. 
This increase cannot be due to (1) in Lemma 4 because has transmitted more bits to the stronger user 
by tu, i.e., Yli=i ^hCi < Yll=i ^liCi- Moreover, this increase cannot be due to (2) in Lemma 4 since rate 
of the weaker user in 5"^ is greater than zero in epoch u, i.e., = h2{Pu, f^u) > ^2(-Pu, rf^) > 0. Hence 
rate of stronger user in cannot increase after epoch s. Finally, rate of the stronger user in 5"^ cannot 
decrease (See Lemma 4) as this would also contradict optimality. Hence, there cannot be two optimal 
schedules with different rate pair vectors. 

As both the power allocation vector and the rate pair vector of an optimal schedule are unique, we 
conclude that the optimal schedule is unique under WUFBC. ■ 
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V. The DuOpt Algorithm 

The problem in [9] which is a special case of Problem 1, where both users' data is available at the 
beginning, was shown to be solved in [9] by the FlowRight algorithm [15]. Along similar lines, we 
develop an algorithm that we call DuOpt for solving Problem 1 in its general form. As a matter of 
fact, DuOpt simply reduces to FlowRight when the given problem instance has all the data arriving at 
t — Similarly to FlowRight, DuOpt starts with any feasible schedule and reduces the transmission 
completion time iteratively. Let the number of epochs and the transmission completion time of the initial 
schedule be A;"^ and T"^ respectively. In each iteration, DuOpt sequentially updates rates and powers of 
two consecutive epochs at a time, i.e., epochs (1, 2), (2, 3), until all epochs are updated. Then, starting 
from the first epoch pair, DuOpt continues with the next iteration. DuOpt stops after N iterations such 
that N = min{n : T""^ -T" < e, i = 1, A;", j = 1, 2}, where T" < T^^p is the transmission completion 
time. A;" < A;"p is the number of epochs used at the end of n*^ iteration and e is a predefined threshold. 

Hereafter, we will briefly outline the local optimizations over epoch pairs. In Theorem 2, it will be 
shown that local optimizations can only improve the schedule. We will also prove that under WUFBC, 
successive iterations strictly improves the schedule unless it is optimal. 

Local Optimizations: Let Ef denote the energy used during the i*'* epoch and 6"^ denote the number 
of bits transmitted to the user during epoch i at the end of n*'* iteration. Suppose that DuOpt is 
at the n*^ iteration and running a local optimization over epoch pair {i.i + 1). The values of V^^ and 
E'^, \/z e {l,2,...,i — 1} have already been found by previous local optimizations. At the end of 
this optimization, E"^ and 6"^ will be determined; E^J^^, E^^2 ^^^d &"(7+i) ^^^^ ^Qstt to new values 
that conserve total energy consumption and data transmission in these epochs. The goal of the local 
optimization is surely to minimize the total transmission completion time of all the packet arrivals. Hence, 
it is logical to minimize the transmission time in the local optimization problem, which results in a gap^ 
if transmission ends before the end of {i + 1)*'* epoch. This gap is used in the next local optimization to 
further reduce the transmission time via transferring bits or energy between epochs + 1) and (i + 2); 
hence, a new gap occurs at the end of the next local optimization. This new gap propagates to the end 
of the transmission resulting a reduction in the total transmission completion time [9]. However, in some 
cases an epoch long gap occurs and this gap is useless for the next local optimization, i.e., energy or data 
transfer between epochs in the next local optimization is impossible because of constraints. In that case, 

gap is a time period with zero power allocation. 
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it is better to just spread the data out till the end of the second epoch in the local problem and minimize 
the energy consumption so that the excess energy can be used to further reduce the transmission time in 
the next local optimization. This leads to two different local optimization functions: time minimization 
and energy minimization. These functions both support the global objective in different ways. Time 
minimization aims to find the minimum amount of time, T^^^^iy to transmit j+i) = b'^^^ + ^"(7+i) bits 
to each user using the energy available in epoch pair + 1), i.e., E^[l^i) = E^~^ + E^+i- On the 
other hand, energy minimization aims to find the minimum energy, E^^^iy to transmit ^'"(jj+i) bits to 
each user in two epoch durations, i.e., + ^j+i, and excess energy, -Bf ~^ + E^~i — EJ^^ ■_^-^y is transferred 
to the (i + 2)*'' epoch in order to conserve energy. Both of the optimizations respect energy and bit 
causalities, i.e., E^ < E{ti) - E^ and b]^ < Bj{ti) - Yfs=\. J ^ 2}. For details of the local 
optimization, see [18]. 

Suppose that all the feasible packets have been transmitted until the end of the i^^ epoch and there are 
still packets to arrive after ti. Then, further minimization of transmission completion time of sequential 
epochs before U will be suboptimal. On the other hand, we can minimize the energy consumption until 
ti and use the excess energy to minimize the transmission completion time. Therefore, utilization of 
energy minimization for local optimizations in Problem 1 is very crucial. If it is guaranteed that current 
schedule uses at least the same amount of energy as optimal schedule until U, DuOpt uses the energy 
minimization function upto i*'* epoch pair and the time minimization function for the rest. In order to 
determine when to switch from energy minimization to time minimization, a Flag is placed at i*^ epoch 
pair. Initially, the Flag is set to zero and DuOpt starts with performing time minimization on epoch pairs. 
During n*'* iteration, if all the feasible bits are transmitted by the i*'* epoch for 3i e {1, 2, /c"}, then 
the Flag is set to i {Flag < i). In the following iterations, energy minimization function is used up to 
i*^ epoch pair. Fig. 2 illustrates the Flag usage and the pseudo-code in Algorithm 1 outlines the DuOpt 
algorithm. 

Theorem 2: Following statements hold: 

1) Successive iterations of DuOpt can only improve the schedule. 

2) DuOpt stops and returns a schedule with {r^,r5^}. 

Proof. 
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Algorithm 1 DuOpt Algorithm 



InitializeO; 

n ^ 0, Flag ^ 0, TO ^ 
repeat 
n++ 

for i = 1 to (A;" - 1) do 

^"i.max ^ E(ti) — X/m=l^^m 

^li.jnax ^ Bi{ti) — X]m=l ^Im 

^2i,max ^ B2{ti) - J2m=l ^2m 
un , un-1 I un-1 
"l ^ "li ^ "l(i+l) 

^2 ^ "2i ^ "2(i+l) 

if i < Flag then 

run un—l un 1 ipn Tpii-—1 TT''"■~^^ — 

L^li' "l(i+l) ' "2i' "2(i+l) ' ' ' J " 

Minimize_Energy {E^K E^-,\ El^^^h^, , b^, e^max ' Ki,max ' f^i,max^ 

13: else 

Minimize_Time {Er\E^-,\b-, 6^ e^^„,, 

end if 

if^H,ma. — &2i,ma. == ^2. Flag < ^ && ^ < /c" - 1 then 

Flag = i 
end if 
end for 

Calculate_T (&T" ) {Calculate current transmission completion time.} 
until T" 



1) Suppose that DuOpt is running its n*" iteration. After the local optimization on i*" epoch pair, 
we obtain {(r",^,),^^} and reset the values of {{r'^^^l^^^r^^^l^^), E]'^,\ El"^^}. If the Flag is not 
placed before i*'* epoch pair, i.e.. Flag > i, then the aim of the local optimization will be energy 
minimization. Following the local optimization on i^^ epoch pair, the excess energy will be transferred 
to E'^~2. In the next local optimization this excess energy is either further transferred or is used to 
reduce the transmit time. On the other hand, if Flag < i, then the aim of the local optimization on 
i*'* epoch pair will be time minimization. After the local optimization the transmission completion 
time of the bits in epochs (i, i + 1) will either be equal to or before the end of the epoch (i + 1). That 
is, a gap may occur within i*^ epoch pair. In the next local optimization, this gap would propagate 
to the (i + epoch [9]. During the n*'* iteration of DuOpt, if a gap occurs or excess energy is 
transferred during local optimizations, then the gap (or the excess energy respectively) will propagate 
to the last epoch pair resulting in an ultimate reduction the transmission completion time at the end 
of the iteration, i.e., T{r^-,r2i) < ^(^S"^' ^S"^)- neither excess energy nor a gap occurs during 
local optimizations, then transmission completion time cannot be decreased and DuOpt will stop by 



11 

definition. 

Both local optimizations are in favor of the next local optimization. Therefore, if in either one of the 
local optimizations a gap occurs or excess energy is transferred then it would propagate till the last 
epoch pair and finally the transmission completion time would decrease at the end of n*'* iteration, 
i.e., T{r'^-,r2i) < T(r"j"^, r""^). If neither excess energy nor gap occurs during local optimizations, 
then transmission completion time would not be decreased and DuOpt would stop. ■ 
2) In Part-1 we have shown that transmission completion time, T{r'^-,r2j), is strictly decreasing in each 
iteration; meanwhile it is bounded below by T^^'^, Therefore, the iterations of DuOpt stop and return 
a schedule {r^, r^}. ■ 

Optimality of the DuOpt algorithm under WUFBC: 

Theorem 3: If WUFBC is guaranteed, the schedule returned by DuOpt is optimal, i.e., T({r^, r^}) — 

Proof. Suppose that DuOpt stopped and returned a schedule {r^,r^} = S^^, with completion time 
^({'^i?) '^2?}) — Let S°^^ be the unique optimal schedule with transmission completion time T°p*. 
We will now prove that S^^ = 5*°^*. Let us suppose S^^ ^ S°^^, then these schedules have to differ in 
either the power allocation or rate allocation (or both). First, suppose — P°'^^,i e {l,2,...,s — 1} 
and 7^ P°P* for some s. We will show that this case is impossible. There are two possible cases for 
epoch s: (i) " > P°p*, (ii) Pf < P°p\ Let us begin with the first case. 

(i) We assumed P^" > P°p*. If P°p* stays constant after epoch s till the end of transmission, this 
would mean that S^^ consumes more energy than 5"°^* until 7"°?*, which would contradict the fact 
that optimal schedule consumes all the harvested energy till the end of transmission. Therefore the 
power of the optimal schedule must increase at the end of epoch (s + m) for some m > before 
the end of transmission. As S^^ has been able to use more energy than the optimal schedule until 
tg+m^ the optimal schedule cannot have run into an energy constraint at tg+m^ hence the rise in the 
power can only be due to a data constraint at ts+m, i-C, all the bits arrived have been transmitted by 
the optimal schedule until tg+m- In order to contradict the assumption that P^" > P°p*, we will now 
analyze the rate assignments for both schedules. First let us focus on the case that both schedules 
use exactly the same rates for the stronger user up to ts+m, i-e-, rf^^ = r°f* , Vi G {s, s + m}. As 
we have shown above, 5'°p* should have transmitted all the bits available until tg+m- However, if we 
compare the weaker user bits transmitted by both schedules until ts+m> we observe that S^^ transmits 
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more bits to the weaker user than 5°?* does, because Ei=r(^?i'' - = YZTiT2t - = 

Ei=r(^2(i^°",rfi") - h2{P^\r'^^))ii > 0. On the other hand, DuOpt respects bit causality, i.e., 
DuOpt does not transmit bits that have not arrived yet, so we reach contradiction. That is, rates 
cannot stay constant up to tg+m, i-e., there is some k e {1, 2, s + m — 1} such that r^^ — r°f* 
for i < A; and rf^" ^ r°f . But we shall now show that this is not possible. First consider the 
case that r^^ < r°^\ From Lemma 4, the stronger user's rate cannot decrease under WUFBC. If 
rf^ — rf^, i e {k, s+m}, then S^^ transmits more bits to the weaker user than 5'°^* does by tg+m-, 

However, at the end of (s + m)*^ epoch, 5"°" cannot send more bits to weaker user because 5"°^* 
should have transmitted all the weaker user bits. Therefore, rf " should increase before tg+m', i-e-> at 
the end of epoch k-\-n, where 0<n< (s + m — A;). We have r^^j^^) ^ ^?(fe+n+i)' hence either one 
of the two conditions in Lemma 4 must hold. Since Y^^i=i{fu^ ~ > 0' until tk+n, -5'°^* has 

transmitted more bits to the stronger user than S^^ does; therefore, all the stronger user's bits arrived 
have not been transmitted by 5'°" at the end of epoch (k + n). Also, r^^^^) ~ ^'^(^k+m''^?{k+n)) ^ 
h2{Pk+ni''^ifk+n)) — 0- Hence neither of the two conditions in Lemma 4 holds and stronger user's 
rate cannot increase at ifc+n, which implies rf^ r°^*. Thus, we are left with the case rf^" > r°^*. 
If r°f* = rll^,i G {k, s + m}, then Yli=^{rf^ — r°f*)^i > , which contradicts the fact that 
DuOpt respects bit feasibility. Hence, stronger user's rate in S°^^ cannot remain constant after epoch 
k. Then we should have r°f* = r°^*, i e {k, k + n} and r°^l_^^^ < r°^l_^^_^^y Since there is an 
increase in the stronger user rate, at least one of the conditions in Lemma 4 should hold at tk+n- 
However, we have T.i=i (PP'' - > and T.i=iiri"-rlf% > 0, which tells us that neither 

one of the conditions in Lemma 4 holds, which implies that this final case is also not possible. 
Hence, we conclude that the case Pf*" > P°p* is not possible, 
(ii) Now consider the case P^" < P°p*. We will prove that this case is also not possible by following 
a similar method to the one in case (i). First, suppose that the power of S^^ increases after s*'^ 
epoch. This increase cannot be due to an energy constraint, since 5"°^* consumes more energy than 
S^^ does until the increase in power. Hence, it should be due to data constraint and under WUFBC 
both user data constraints should be active. That is, S^^ transmits all the feasible data until the 
increase in power. This implies that while consuming less energy, S^^ transmits at least the same 
number of bits than 5"°^* does, which contradicts the optimality of 5"°^*. Thus, power of S^^ cannot 
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increase after epoch s. Also, it cannot decrease in time, otherwise a local optimization results in 
either a gap or excess energy that propagates till the end of the schedule and transmission duration 
decreases. Therefore, we power of S^" should stay constant after epoch s until T°p*. Now, we will 
analyze the rate assignments for both schedules. Let the transmission of 5"°^* end in epoch (s + m) 
for m > and suppose that rf^ — r°f* Vi < k, < k < (s + m). At the A;*'* epoch there 
are three possible cases: rf^" > r°^\ rf^" < r°^* and rfj^ = rlf". We will first consider the case 
rj^y > rll^ and prove that this is not possible. Let rf^ > r^^* and consider the rate of the stronger 
user in 5"°^* after k*^ epoch. It cannot stay constant until T°^^, because it contradicts the fact that 
S°P^ transmits all the feasible bits before T°p*, i.e., 'ZZTi^ii' - rT)^i > 0- Since the stronger 
user's rate in 5"°^*^ cannot decrease by Lemma 4, it should increase at the end of epoch {k + n) for 
0<n<{s + m-k), i.e., rl^^^^ < r^^l^^^.y However, we have E?=i"(^h" - rT)^i > and 
^2T+n) = HPkllr'lfUn)) > h2{Pk+n.r%+n)) > which implics that none of the conditions in 
Lemma 4 holds and the stronger user's rate in S^^ cannot increase after epoch k. Hence, we conclude 
that rf^ i> . Now we consider the case rf^" < r^^\ Suppose that the stronger user's rate in 5"°^* 
increase at epoch (/c + n) for < n < (s + m — /c). This increase in stronger user's rate requires that 
at least one of the conditions in Lemma 4 should hold. However, we have X]f=r(^B" ~ > 

and r^ffc+n) = ^zlPfcTn'^Ufc+n)) > ^2(Pfc+"n, r°(^+„)) > 0, SO the stronger user's rate in S^"" cannot 
increase, i.e., rf^" r^^*. Since the stronger user's rate in S°^^ cannot decrease by Lemma 4, it 
should stay constant until T°p*. 

Thus far we have shown that if 5"^" is different than 5"°^*, then 5"^" cannot have higher power level 

than 5'°P'^ until T°p*. Moreover, if power level of becomes lower than that of 5°^*, then it should 

stay constant until T°p* and if the stronger user's rate of 5"°" becomes lower than that of 5"°^*, then 

it should stay constant until ^"p*. These results are shown in the general case in Fig. 3. 

Now, let 62 " and be the number of bits transmitted to the weaker user till ts+^+^-i by S^^ and 
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till T°P* by S'°P\ respectively. Then, we have 

s-\-m-\-l s-\-m-\-l 

i=l i=l 

S—1 / s-\-m-\-l \ s+m 

«— \ i=5 / i=s 
/ s+m+i \ s+m 



l—S 



/ s+m+i \ s+m ^ s+m 

z=s /-^i=s i=s 



/ 5+m+/ \ / s+m ^ \ s+m 

> ( E ^0 ( E ^..lTr^/^2(i^°^\r-) - E ^MPF^,rT) 

L^i—s / i=s 



\ l — S / \ l — S ^-^l — S 

s+m s+m 



= E ^MP^^\ r?/) - E ^MP?^\ rlT) 

i—s i—s 

s+m 

= E^*(/^2(^°'*'^?")-^2(/^°''*,rr) 

>0 

> (5) 

From (5), S^^ transmits more bits to the weaker user than 5"°^* does, then this final case also cannot 
happen. Therefore, we conclude that the schedule returned by DuOpt cannot be different than the 
unique optimal schedule, i.e., 5*°" = 5*°^*. ■ 



VI. Numerical Example 

Consider a two-user AWGN BC with IKhz bandwidth and A^o = 10"^^ Watts/Hz. Path loss factors 
on the links of stronger and weaker user are assumed to be 70dB and 75dB, respectively. Amounts and 
instants of energy harvests and bit arrivals are depicted in Fig. 4. Under these circumstances DuOpt 
algorithm is run and final schedule is calculated as drawn in Fig. 4 

VII. Conclusion 

With the new advances in energy harvesting technologies, optimization of communication systems that 
depend on renewable energy resources has emerged as an important problem. A large body of recent 
research effort in the field has focused on transmission scheduUng policies [9], [10], [6], [7], [12], [13]. 
This paper continued the work on to the solution of the previous formulation of the problem stated in [9] 
and [10]. In particular, it aimed to find the power and rate allocation policy in a Broadcast Channel with 
two users, that minimizes the total transmission completion time of data that becomes available at arbitrary 
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points in time, with energy that is harvested at arbitrary points in time. It should be noted that this Une of 
work relies on offline problem formulations where exact knowledge of data and energy arrival events is 
assumed. The approach is useful for obtaining structures and guidelines as well as benchmarks for online 
problems which may be more applicable in practical transmision scenarios. 

In this study, structural properties about the solution of the general optimal offline broadcast packet 
scheduling policy for an energy harvesting broadcast channel have been established. The uniqueness of 
the optimal policy under the weak-user-full-buffer condition (WUFBC) has been shown. An iterative 
algorithm, DuOpt, that returns a feasible schedule which possesses the same structural policies that the 
optimal was shown to have is devised. Moreover, DuOpt was shown to obtain the unique optimal policy 
under WUFBC. 

Proving the uniqueness of the optimal poUcy and optimality of DuOpt in general are two goals of 
ongoing work. The direction set for future work related to the problem presented in this paper also 
includes online problems. 

VIII. Appendix 

A. Proof of Lemma 2 

The claim is that power is non-decreasing with epoch number i. Equivalently power is non-decreasing 
in time. To show this, we will argue that given a schedule in which power decreases at some time ti, 
this schedule can only be improved by equating the power levels before and after ti. Consider a time 
interval (ri,r2), so that power is constant at Pi > during {ti^U), and at P2 < Pi during (ti,r2). As 
illustrated in Fig 5, let t — T2 — ti, and the lengths of the constant-power slots be /3t and (1 — /3)t. Denote 
the rate pairs in the 1** and 2"'^ slots as (rn, r2i) and (r^, r22), respectively. We will show that keeping 
the total consumed energy constant, and transferring some amount of energy AE from the first slot to 
the second such that power levels are reallocated closer to each other, the sender can transmit at least 
the same number of bits within the same duration. Let us denote the average rate of the weaker user as 
r2 = l3r2i + (1 — /9)^22- Provided r2i > 0, the sender could transfer some energy and some of user 2's 
bits from the first epoch to the second while keeping user I's rates rn and constant. As energy and 
bits are simply being deferred for later use, this operation does not violate feasibility. Specifically, let 

P; = Pi - (1 - P)AP , P2 = P2 + PAP (6) 
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such that the new power allocation to the slots is (P{, P2) satisfying P2 < P2 ^ Pi ^ Pi- With this new 
allocation, the weaker user's rate in the first slot is h2{Pi,ru) < h2{Pi,rn) = ^21 > 0. Its new new 
average rate over the duration of t is: 

r=2 = /i2(P;,rn)/3 + /i2(P2,^i2)(l-/3) 

> /i2(Pi,rii)^ + /i2(P2,ri2)(l-^) (7) 

= ^2 

This is shown by straightforward application of the properties listed in section II. 

In the remaining case which is r2i = 0, we know that rn > must hold (as Pi > 0.) In this 
case, the allocation can similarly be improved by bringing power levels closer and transferring some of 
the first user's bits to the right, while keeping the rate allocation of the weaker user unchanged. Let 
n — Prii + (1 — /5)ri2 be the average rate of the stronger user over the duration t. After the reallocation, 
the average rate of the stronger user becomes 

r=i - /ii(P;,r2i)/3 + /ii(P2,r22)(l-/3) 

> hi{Pi,r2i)/3 + hi{P2,r22){l-/3) (8) 
= fi 

(8) follows from the fact that 

hi{P'i,r2i)/3 + /ii(P2,r22)(l - P) - hi{Pi,r2i)/3 - /ii(P2, r22)(l - /?) > (9) 
for all /3 = {0, 1} with equality achieved at /3 — 0,1. 

q{P) = /ii(Pi-(l-/3)AP,r2i)/3 + /ii(P2 + /3AP,r22)(l-/3) 

-hi{Pi,r2i)p - /ii(P2,r22)(l - P). (10) 

We can show that (9) holds by proving q(/3) is concave in /3. 

The 1** and 2"*^ order derivatives of q with respect to /3 are the following^ 



^hi^ and hi^^ represent the first and second order partial derivatives of hi with respect to P, respectively. 
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dq 
dp 



9^2 



= /ii(Pi - (1 - /3)AP,r2i) +^{/iijPi - (1 - /3)AP,r2i)(AP)} 
-/ii(P2 + /5AP, r22) + (1 - /5) {/ii. (P2 + /3AP, r22)(AP)} 
-h{Pi,r2i) + h2{P2,r22) (11) 

= 2(/iijPi - (1 - /3)AP,r2i)(AP) - /iijP2 + /3AP,r22)(AP)) 

" V ' 

<0 

+/3Ui..(A-(l-/3)AP,r2i)(AP)^ 



<o 



+(l-/9)Ui..(^2 + ^AP,r22)(AP) 



2 



<0 



< (12) 

According to the properties listed in section n, (12) always holds if r2i > r22. Hence q is concave in 
, if 1^21 — 0. We conclude that a policy that contains a drop in power level is sub-optimal. ■ 



B. Proof of Lemma 3 

To reach contradiction, suppose that power increases at time U (Pj+i > Pi). We will show that if none 
of the conditions (a), (b) or (c) hold, then it is possible to improve the schedule by transferring some 
energy from the (i + 1)*'* epoch to the i*'^. Assuming the i*^ epoch length is /3t and the + is {l—/3)t, 
after bringing power levels closer, 

P-^Pi + PAP , P^+i = P,+i - (1 - ^) AP (13) 

Observe that if we treat Pi as P2 and P^+i as Pi, then (13) becomes identical with (6). This implies 
that at least the same number of bits could be transmitted to weaker user, if we can bring power levels 
closer while keeping the stronger user's rates constant. In addition to this, the allocation could also be 
improved by bringing power levels while keeping the weaker user's rates constant, in case r2{i) > ^2(1+1). 
Consequently, equations (7) and (8) hold. 

It is straightforward that we cannot bring power levels any closer when condition (a) holds, due to the 
energy causality constraint. Secondly, it also doesn't yield a better schedule, if we can not transfer data 
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from the latter epoch to the former(condition (b)). When it is possible to transfer some positive amount 
of energy from the {i + epoch to the i*'^, as shown in (7), we can always improve allocation while 
keeping the rates of the stronger user the same. Although bringing power levels closer while keeping the 
rates of the stronger user the same is not feasible in case weaker user's data constraint is active, we may 
still improve allocation as proved in (8). Nevertheless, this time we require r2(i) > r2(j+i). As rate of the 
weaker user can only rise upon a data arrival for the weaker user in case weaker user's data constraint is 
active, condition (c) describes the last case that we may not improve allocation by bringing power levels 
closer. 

We have thus shown that this set of three conditions contains all the cases in which power can rise, if 
none of these hold, then power cannot rise. It is straightforward to show that this set cannot be further 
reduced by finding counterexamples for the claim that if any one of (a), (b) or (c) is satisfied, then power 
may rise at time ti. ■ 

C. Proof of Corollary 

1) Suppose that power increases upon a bit arrival for the weaker user occurring at ti. As there is no 
energy constraint at ti, bringing power levels closer does not contradict with the energy causality in 
this case. This implies that conditions (b) or (c) stated in Lemma 3 must hold. However, we know 
that there is a data arrival for the weaker user at ti, so if (b) were true, then (c) would be true as 
well. Therefore, condition (c) holds in either case. ■ 

2) Suppose that power increases upon a bit arrival for the first user. With similar reasoning to part- 
1, condition (a) of Lemma 3 cannot hold. As there is also no data arrival for the weaker user, 
condition (b) must be satisfied. ■ 

3) As there is no data arrival at the time when power increases, the only possibility that power increases 
upon an energy harvest is condition (a) of Lemma 3. ■ 

D. Proof of Lemma 4 

Suppose that ru < rm^iy One can find a better schedule by bringing the rates of the stronger user 
closer by Lemma 6 of [9]"^. Therefore, the stronger user's rate cannot decrease. However, the stronger 
user's rate may increase because it may be against to either bit or energy causality to transfer some 

''Lemma 6 of [9] originates from the observation in [10] (See Lemma 4 and Corollary 1) that there is a cut-off level for the total power, 
below which no power is assigned to the weaker user. 
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Stronger user bits from epoch i + 1 to i. Firstly, it is against bit causality to transfer some stronger 
user bits from epoch i + 1 to epoch i, if the first condition holds. Secondly, if the second condition 
is satisfied we cannot bring stronger users rates closer to each other as it would violate energy causality. ■ 

E. Details of the Local Optimization 

Consider the local optimization problem given in Fig. 6, where 5^ is the data arrival for the i*'* user, 
Ej is the energy harvest at the beginning of j*'* epoch for i, j e {1, 2}. Ti is the length of the first epoch 
and T2 is the transmit duration in the second epoch. Let Pij and Vij be the power and rate assigned to 
the i*^ user during j*'' epoch after optimization. The energy and data causality constraints for the local 
optimization problem are as follows: 

El > Pn-Ti 
Bn > rn-Ti 
B21 > rsi-Ti. 

The structure of the solution changes if either one of the constraints satisfied with equality. Since there 
are 3 different constraints, after optimization one will encounter one of the 2^ = 8 results. We have studied 
all 8 cases and derived solutions to each one of them for both energy minimization and time minimization 
functions. In each case, the solution can be calculated analytically for energy minimization or it can be 
found iteratively for time minimization functions. Before starting the optimization, if one already knows 
which constraints should be satisfied with equality, then the result could be obtained solving just that 
case. Otherwise, one should compute the results for each case and then select the best one^ that respects 
energy and bit causalities. Next, we will analyze one of the cases and present the algorithms for both 
energy and time minimization functions. The analyzes of the remaining cases can be done in a similar 
fashion [18]. 

Local Optimization when only stronger user data constraint is active: One of the possible structures 
of the local optimal solution is when only the stronger user data constraint is active. In this case, total 
transmit power level in local optimal solution should be constant (See Lemma 3). Under this condition, 
a possible illustration of optimal power allocation is depicted in Fig. 7. Since the total power should stay 

^Best result is the one that consumes minimum energy for the energy minimization function and one that has the minimum transmit time 
for the time minimization function. 
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constant during transmission, we should have P11 + P21 = ^2 + ^22- Moreover, we should have Pu < P12 
since only the stronger user bit causality is met. 



For a given transmission completion time Ti + T2, the optimum schedule that minimize energy 
consumption can be found as follows: 

Pii = -(2^1 -1) , Pi2 = — (2^2 -1) and P22 ^ Pii + P21 - P12 (14) 

Sl Si 
B2 = ^21 + 522=-l0g,(l + ^^^^-^) + -l0g,(l+^^^^-^). (15) 

From (14) and (15) we derive 

P21 = -Pii + - [(2'^H^ii«2 + (jY'{Pi2S2 + cjY') ^ - a'l . (16) 

S2 L J 

Then, total energy consumed in two epochs is calculated by 

^mm = (Al + P2l)(ri+T2). (17) 



Minimum energy to transmit Bi = Bu + B12 and B2 = -B21 + ^22 bits to the users in two epochs, 
Emin, can be calculated by setting T2 as the length of the second epoch in (17). 



Substituting (14) into (15) and arranging terms we obtain 



B. 




Ti + T2. fE^ + E2 , (72\ 



2 °" V ^1 + 7^2 S2J 



The first and second order derivatives of B2 with respect to T2 are as follows: 
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dB2 1 , fEi + E2 



T. + n sj 21n(2)(l + g|i±g) 



-- log2 -2 ^2 + a' + 

2 \si S1S2 J I2 I I si-s2 2~^ 

2S, 



92^2 1 21n(2)i?22^2 -2 



2in(2)(r, + r2)(i + gg±f)' 73(1 + ^2-^)^ 



(18) 



< 



As shown in (18), second derivative of B2 is always negative for si > S2, which implies that B2 is a 
strictly concave function of T2. As T2 goes to infinity, B2 is as follows. 



hm B2 = — — log2 —(2 ^1 - 1) + 1 B12 + 



T2-700 2 ^^\si^ ' J Si 21n(2)(72 

> —00 

Since B2 is a strictly concave function of T2 and goes to a finite number as T2 goes to infinity, B2 is 
an increasing strictly concave function of T2 and there is a unique B2 for each value of T2. 

In time minimization we have Ei, E2, Ti, B2 = B21 + -B22> -B21, -B22; Si and S2 are constant terms and 
T2 is transmission time within the second epoch which is to be minimized. Since B2 is an increasing 
concave function of T2, we iteratively find T2 that sends exactly B2 bits to weaker user by using bisection 

method. Minimum time to transmit Bi = Bu + B12 and B2 bits to the users in these two epochs can be 

(2) 

calculated by — T1 + T2. Algorithm 2 presents a pseudo-code of time minimization algorithm for 
the case that only the stronger user bit causality event occurs. 
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Fig. 1. Example for (a) a sequence of energy and data arrivals, (b)-(c) the corresponding E(t),B{t), (d) the schedule P(t) and {ru, r2i}- 
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Fig. 2. Illustration of the Flag and local optimizations, where all the feasible bits have been transmitted until the end of 5*'' epoch; hence, 
a Flag is set to 4"* epoch pair, i.e., (4,5). Energy minimization is performed upto epoch pair with the Flag and time minimization is 
performed for the rest. 
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Fig. 3. Illustration of the final case in the proof of Theorem 3. 
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Fig. 4. (a) An illustration of energy harvest and bit arrival sequences, (b) Final schedule calculated by DuOpt, where Ps is the power used 
for transmission to stronger user and Pt is the total power used by the transmitter. 
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Fig. 5. Illustration of the transmission scheme used in proof of Lemma 2. 
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Fig. 6. Illustration of the local optimization problem with two epochs. Bij represents the data arrival for the j user at the beginning of 
the i*** epoch. Similarly, Ei represents the energy harvest at the beginning of the i*'' epoch. 
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Fig. 7. Optimal power allocation for the stronger and the weaker user if only the stronger user data causality constraint is active. Dark 
shaded levels represent the stronger user's power levels, whereas the light shaded ones represent the weaker user's power levels. The total 
transmit power stays constant. 



